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^ ■ ABSTRACT 
O 

^ I We report on the discovery of a faint (My ~ —10.6 ± 0.2) dwarf spheroidal 

galaxy on deep F606W and F814W Hubble Space Telescope images of a Virgo 
, intracluster field. The galaxy is easily resolved in our images, as our color mag- 

O ! nitude diagram (CMD) extends > 1 magnitude beyond the tip of the red giant 

branch (RGB). Thus, it is the deepest CMD for a small dwarf galaxy inside 
a cluster environment. Using the colors of the RGB stars, we derive a metal 
\Q I abundance for the dwarf of [M/H]= —2.3 ± 0.3, and show that the metallicity 

^ I dispersion is less than 0.6 dex at 95% confidence. We also use the galaxy's lack 



O ' -"^ Department of Physics & Astronomy, Youngstown State University, Youngstown, OH 44555; prdur- 

-fi . rell@ysu.edu; jjfeldmeier@ysu.edu 

' ^Department of Astronomy & Astrophysics, Pennsylvania State University, University Park, PA 16802; 

bwilliams@astro.psu.edu; rbc@astro.psu.edu; steinn@astro.psu.edu 

' '^The University of Texas, Department of Astronomy, 1 University Station C1400, Austin, Texas 78712; 

■ ted@astro.as.utexas.edu 

■^WIYN Observatory, 950 North Cherry Avenue, P.O. Box 26732, Tucson, AZ 85726; jacoby@wiyn.org 

^Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218; ferguson@stsci.edu 

^Department of Physics & Astronomy, University of Leicester, Leicester LEI 7RH, UK; nrt3@star.le.ac.uk 

^European Southern Observatory, Karl-Schwarzchild-Str. 2, 85748 Garching, Germany, marnabol@eso.org 

^Max-Planck-Institut fuer Extraterrestrische Physik, P.O. Box 1312, D-85741 Garching, Germany, ger- 
hard@exgal . mpe . mpg.de 

^Instituto de Astrofisica de Canarias, C/ Via Lactea, s/n, 38205, La Laguna, Tenerife, Spain; jal- 
fonso@ll.iac.es 

^°Mount Stromlo Observatory, Research School of Astronomy and Astrophysics, Mount Stromlo Observa- 
tory, The Australian National University, AGT 0200 Australia; kcf@mso.anu.edu.au 

^^NSF Astronomy & Astrophysics Postdoctoral Fellow 



-2- 



of AGB stars and the absence of objects brighter than M^oi ~ —4.1 ± 0.2 to 
show that the system is old (t > 10 Gyr). Finally, we derive the object's struc- 
tural parameters, and show that the galaxy displays no obvious evidence of tidal 
threshing. Since the tip of the red giant branch distance ((m — M)o = 31.23±0.17 
oi D = 17.6 ± 1.4 Mpc) puts the galaxy near the core of the Virgo cluster, one 
might expect the object to have undergone some tidal processing. Yet the chem- 
ical and morphological similarity between the dwarf and the dSph galaxies of the 
Local and M81 Group demonstrates that the object is indeed pristine, and not 
the shredded remains of a much larger galaxy. We discuss the possible origins of 
this galaxy, and suggest that it is just now falling into Virgo for the first time. 

Subject headings: galaxies: dwarf — galaxies: stellar content — galaxies: abun- 
dances — galaxies: clusters: individual (Virgo) 



1. Introduction 



An understanding of the processes that affect the formation, evolution, and sometimes 
destruction of dwarf galaxies is critical to our overall picture of galaxy formation. Since dwarf 
galaxies are the most common type of galaxy in the Universe, any model of galaxy formation 
is incomplete if we cannot understand the se objects. Moreover, under the current paradigm 



of hi erarchical structure formation (e.g., iKauffmann. White fc Guideronil Il993l : ICole et al. 



20001 ) . massive galaxies are formed from the mergers of less massive objects. Thus, dwarf 
galaxies are an important building block of the universe, and an understanding of their 
properties will help illuminate the processes of galaxy formation more generally. 

The faintest examples of d warf galaxies are the dwarf spheroidals (dSphs). Using the 



nomenclature of iGrebel et al.l (120031 ). these are low- luminosity (My > —14), low central 
surface brightness (/xy > 22 mag arcsec"^) objects with shallow radial profiles and stellar 
populations that are dominated by old and intermediate-aged star^. Most of what we know 
about dSphs comes from observations of ~ 20 galaxies in the low-density environment of 
the Local Group. From the data, it appears that dSphs are fairly homogeneous in nature, 
with small, but dense dark niatter halos and low metallicities that are in proportion to their 
luminosity (for a review see iMated Il998l ) . 



Much less is known about dSphs in the cluster environment. Surveys in Virgo and For- 



^dSph galaxies are the faint end of the dwarf eUiptical [dE] galaxy sequence, and when we note dE galaxies 
we are including those gas-poor dwarf galaxies more luminous than My = — 14 
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nax have successfully found large numbers of dE and bright dSph galaxies (|Bingge 



19851 : iFergusonlllQSa IPhiUipps et al.lll998l : IXrentham fc Hodgkinll2002l : ISabatini et al. 



i et al. 



20031 ) 



but there is still considerable u ncertainty in just how many of these lowest-luminosity objects 
exist 



le.g. 



Hilker et al.ll2003l ). The scarce information we do have on cluster dSph galax- 



ies suggests that the ir properties are similar to those of Local Group dwarfs (IHilker et al. 



20051 : ICaldwellll2006l ). If correct, then such a result is somewhat surprising, since the tidal 
forces of the clust er environment are expected to take their toll on these objects. Indeed, 
Bekki et al.l (120011 ) have shown that dSphs can be stripped of their outer stars in only a few 
Gyr. This would imply that both the morphology and chemistry of cluster dSphs should 
differ significantly from their Local Group counterparts. 

Until recently, abundance studies of cluster dSph galaxies (based on photome try of 
individual stars) have been impractical. For example, although [Harris et al.l (119981 ) were 
able to resolve the brightest red giant branch (RGB) stars in a Virgo cluster dE,N galaxy 
using a 32 ks exposure with the Hubble Space Telescope's WFPC2 camera, these data were 
limited to a single filter (F814W). However, with the advent of the Advanced Camera for 
Surveys (ACS), it is now possib le to measure the broad-band colors of red giants in these 
systems. Indeed, ICaldwelll (120061 ) has recently used the ACS to resolve the stellar populations 
of four dE/dSph galaxies in the Virgo cluster, and obtain color-magnitude diagrams of their 
stars down to J ~ 28, making it the deepest study of individual stars in cluster dwarf 
galaxies. Their data appear to confirm that the dE/dSph galaxies of Virgo obey the same 
luminosity-metallicity relation as the dE/dSph galaxies of the Local Group. 

Here, we present an extremely deep F606W and F814W study of a faint {My ~ 
— 10.6) dSph that was serendipitously found during a survey of intracluster stars in Virgo 
( IWilliams et al.l l2006al ) . We use the data, which reach I ~ 28.4, to show that the galaxy 
is remarkably similar to undisturbed dSphs in the Local Group, both in morphology, and 
in stellar content. In Section 2, we describe both our data reductions and our artificial star 
experiments, and present a color-magnitude diagram of the dSph's red giant branch stars. In 
Section 3, we use these data to derive the system's distance as well as its stellar population. 
We show that the galaxy has a mean metallicity that is very low ([M/H] ~ —2.3 ± 0.3), 
and is composed entirely of stars older than 8 to 10 Gyr. We also derive the galaxy's struc- 
tural parameters, and show that its central surface brightness and core radius are typical 
of dSphs in the Local Group. In Section 4, we compare the galaxy's properties to those of 
Local Group objects, and attempt to investigate the effects of environment on its history. 
We show that the galaxy displays no obvious evidence of tidal disruption, and has a mean 
metallicity appropriate for its luminosity. We conclude by considering the possible origin of 
this object. 
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2. Observations and Data Analysis 



Between 30 May 2005 and 7 June 2005 we used the Advanced Camera for Surveys on the 
Hubble Space Telescope to obtain deep F606W and F814W images of a Virgo intracluster 
field (a(2000) =12:28:10.80, 5(2000) =12:33:20.0, orientation 112.58 degrees). This field, 
whic h is projected ~ 180 kpc from M87 (assuming a mean Virgo distance d = 16 Mpc; 
e.g., iJerjen et al.ll2004l ) was chosen to be near the cluster center, but away from any known 
object: the closest cataloged galaxy is the dwarf elliptical VCC 1051, which is ~ 5' (24 kpc) 
to the northwest; the nearest large galaxy, M86, is over 170 kpc away. 

Our F814W (/-band) data consisted of 22 exposures, with 26880 s of integration time; 
the F606W (wide V-band) observation included 52 exposures totaling 63440 s. These data 
were co-added, and re-sampled using the multidrizzlcltask within PyRAI^ (jKoekemoer et al 
2OO2I ) to produce images with 0'.'03 pixel"^. The details of these reductions, and a n image of 
the field illustrating its position in the cluster is given by Williams et al.l (j2006a ). 



Immediately after data acquisition, an inspection of our images revealed the pres- 
ence of a previously unknown galaxy in the field. The object, whose center is located at 
a(2000) =12:28:15.5, 5(2000) =+12:33:37.0 (uncertainty ~ 0'.'2) is ~ 15" in extent and 
clearly resolved into stars. A color image of the object is shown in Figure [H 

Our point source photometry was similar (but not identical to) that performed by 
Williams et al.l (l2006al ). To avoid dealing with variations in the ACS point spread function 
(PSF), we began by limiting our photometric reductions to a 2200 x 1600 pixel (66" x 48") 
region surrounding the galaxy. We then chose three bright, isolated stars on our F606W 
image, and four bright isolated stars on the F814W image t o define the P S Fs, a nd used 
DAOPHOT II and ALLSTAR to perform the photometry JstetsonI Il987l . Il992h . Two 
DAOPHOT II/ALLSTAR passes were performed to detect as many of the stars as possible. 
The F606W and F814W datasets were then merged with a 1.5 pixel matching criterion to 
create a preliminary object catalog of ~ 1000 objects. Finally, we used the DAOPHOT x 
and r_2 image parameters to rern ove background galaxies and severely blended stars from 
the list JstetsonI Il987l : lKronlll98ol l 



Our instrumental magnitudes we re placed on the VE GAMAG system using the pre- 
scription and zero points provided by ISirianni et al.l (120051 ). To obtain the offset between 
our ALLSTAR magnitudes and the VEGAMAG system, we used the profiles of the PSF 



^multidrizzle is a product of the Space Telescope Science Institute, which is operated by AURA for NASA, 
http: / / stsdas. stsci . edu/ pydrizzle / multidrizzle 

■^PyRAF is a product of the Space Telescope Science Institute, which is operated by AURA for NASA. 
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stars on each image; the rms scatter in this cahbration was ~ 0.02 mag. We note that since 
the region under consideration is small (only ~ 8% of the entire ACS field), differential 
effects associated with charge transfer efficiency were negligible (~ 0.01 — 0.02 mag), and 
were ignored. 



2.1. Artificial Star Experiments 

To ascertain the photometric uncertainties and incompleteness in our data, we used the 
traditional method of adding simulated stars of known brightness to the science frames and 
re-reducing them following the exact same procedures as for the original data. To each frame, 
we added 9000 stars (300 runs with 30 stars added per run) to an elliptical region centered 
on the dwarf galaxy, and re-executed our two-pass DAOPHOT II/ALLSTAR photometric 
procedure, including the rejection of non-stellar sources, and the merging of the datasets. By 
limiting the number of fake stars to 30 per run, we did not significantly alter the crowding 
conditions on the images; by defining the mean F606W— F814W color of our artificial stars 
to be roughly equal to that of the observed objects (F606W— F814W = 1.0), we ensured 
that our experiments adequately reproduced losses during the catalog merging and image 
classification processes. Finally, by defining the stars' luminosity function to be a rising 
exponential between F814W = 22 and F814W = 29, we improved our statistics at the faint 
end of the luminosity function, where the uncertainties are largest. The result of these 
experiments were the photometric completeness function f{m), defined as the ratio of the 
stars recovered to stars added regardless of the magnitude that was actually measured, A(m), 
the mean shift in magnitude of the stars (A(m) = rriout—^in) , and cT(m), the dispersion about 
this mean. All of these functions were computed in 0.5 mag bins, to ease the computational 
requirements. 

To ascertain each filter's limiting magnitude (defined as the 50% completeness fraction), 
we performed similar artificial star tests using a variety of F606W— F814W input colors. 
The results from these experiments are illustrated in Figure [21 From these experiments, we 
obtained rrium = 29.1 ±0.1 in F606W and 28.4 ±0.1 in F814W, where the errors refiect both 
the rms scatter due to the use of different colors, and the uncertainty associated with varying 
levels of c rowding within t he gala xy. As expected, these limits are slightly brighter than those 



found by IWilliams et al.l (j2006al ). whose analysis focused on the uncrowded regions of the 
frame. Moreover, our classification criteria kept only those objects with the best photometric 
accuracy: obvious galaxies and extremely blended objects were removed from the catalog. 
As result, even at F814W = 27, ~ 10% of the artificial stars were not recovered. Fortunately, 
since we are not making any inferences based solely on the absolute number of stars detected. 
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these missing objects do not affect our analyses in any significant way. 

A color-magnitude diagram for our artificial stars is shown in Figure O The figure 
illustrates that the photometric errors increase with magnitude, so that by mum, the un- 
certainties are ~ 0.2 mag. The figure also shows that there is a small blueward shift at 
faintest levels; this is due to the uneven depth of the two images. (At the faintest levels, 
our F814W photometry only detects objects whose photometric errors are in the positive 
direction.) Fortunately, down to the limiting magnitudes of the frames, this color shift is 
< 0.1 mag. These corrections, while small, will be accounted for in the following sections. 



2.2. Color Magnitude Diagram 



The color-magnitude diagram for all the stellar objects in our field is displayed in the 
left-hand panel of Figure |H The diagram shows that most objects lie in a sequence near 
F606W— F814W ~ 1 that extends up to F814W ~ 27.1. However, there is a significant 
amount of scatter about this line. This scatter comes from two sources: the RGB stars that 



pervad e Virgo's intracluster space (iFerguson et al.lll998l : iDurrell et al.ll2002l : IWilliams et al. 
2006al ). and, to a lesser extent, unresolved background galaxies. To minimize the effects 
of these contaminants, we identified a subsample of stars within the galaxy's F814W ~ 
26.5 mag arcsec"^ isophote. This subsample, which is defined via an ellipse with semi-major 
axis 8'.'1, an axis ratio of b/a = 0.5, and a position angle of 43° (see Section 3.4) should have 
minimal contamination: based on star counts in the rest of our field, only ~ 16 out of the 
171 objects brighter than FSIAW = 28.4 should be contaminants. Thus, we can use the 
data from this 103 arcsec^ region to form a "clean" color-magnitude diagram of the galaxy's 
stars. The elliptical region which defines our subsample is shown in Figure [H the right-hand 
panel of Figure H] shows its CMD. 



3. Analysis 

The CMD in Figure H] displays the hallmarks of a red-giant branch population with 
F814W > 27. There are few stars brighter than this limit, and no evidence of any bright 
blue stars arising from a young stellar population. Since the data clearly extend more than 
1 mag down the RGB, we can use the diagram to investigate both the distance and metallicity 
of the galaxy. 
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3.1. TRGB Distance 



The tip of tlie red giant brancli (TRGB) in the /-band has repeatedly been shown to 
be an excellent distance indicator for old, metal-poor populations. In stellar sys tems with 



[M/Hj < —1, the TRGB is essentially independent of both age and metallicitv fe.g.jLee et al. 



1993; 


Sakai et al 


1996; 


Ferrarese et a 


2000; 


Bellazzini et al. 


2001 


Mouchine et al. 


2OO5I; 


Caldwell 


2006) 


, so all that is required is 



2004 



end truncation of t he RGB and an estimate of foreground reddening {E{B — V) = 0.025; 



Schlegel et al.l 119981 ). Our only limitation is the small size of the galaxy, which restricts the 



number of stars available for analysis. 

Although foreground/background contamination in our 'dwarf-only' CMD (right panel 
of Figure Hj) is small, to make our TRGB detection as unambiguous as possible we have 
further restricted our sample to objects with F606W— F814W < 1.3. Stars redder than this 
are not on the metal-poor RGB and are unlikely to be members of the galaxy. This makes 
very little difference to the analysis, since it excludes only 3 objects from consideration. 
Moreover, since the RGB tip for red, metal-rich stars is fainter than that for blue objects, 
eliminating these stars should not affect our distance determination. We note that the 
ACS F814W fi lter b andpass is very similar to that of the traditional /-band: according to 
Sirianni et al.l (120051 ). the transformation between these two systems is / = FSIAW — 0.008 
(for stars with F606W-F814W=1) in the VEGAMAG system. 

Visual inspection of the Figure H] shows a rather sharp transition at F814W ~ / ~ 27.1, 
albeit with a small number of stars. While the increased Poissonian noise from small numbers 



can lead to a possible systematic bias in /-r^f^r^R. lMadore &: FreedmanI (119951 ) have shown that 
with at least ~ 100 stars in the top magnitude of the RGB, such biases are small. Since our 
analysis includes ~ 130 objects, we are safely above this threshold. 

Another possible source of systematic error is image crowding: in dense fields, stellar 
blends can make the RGB tip appear brighter than it should be. However, in our case, this 
effect should not be important. Our artificial star experiments demonstrate that the mean 
magnitude shift at F814W ~ 27 is very small, less than 0.01 mag. This result, which is 
confirmed by additional artificial star experiments with F814W < 27.1, is due principally to 
our image classification criteria, which rejects the most heavily blended objects (i.e., those 
most likely to bias the results towards brighter magnitudes). 

A quantitativ e estimate of the RGB tip location can be derived using the Sobel edge- 
detection filter of iLee et al.l (119931 ). This technique, which was originally used on binned 
functions, has since been modified to work on discrete stellar samples by treating each 
object as a Gaussian distribution with a dispersion equal to the expected photometric error 
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( ISakai et al.l Il996l ). We note that more rigorous tech niques have been eraployed to derive 



IxBnBi such as the method of maximum-hkehhood (IMendez et al 



2002 



Mouchine et al. 



20051 ) . and the second derivative of the luminosity function (jCioni et al.ll2000l ). However, 
maximum-likelihood methods are better for datasets where the RGB power-law slope is 
observed for over a magnitude. In our case, photometric incompleteness and image crowding 
make the use of the technique problematic. Moreover, the small number of stars limits the 
reliability of second derivative calculations. Because our background contamination is so 
small, the simpler Sobel edge-detection filter should be adequate for our purpose. 

To employ the Sobel-edge detection method, we treated each star as a Gaussian distri- 
bution with a mean equal to its observed magnitude (plus the small Am shift predicted by 
our artificial star experiments), and a dispersion, (^(m), equal to the photometric uncertainty 
at that magnitude. We then co-added the Gaussians, to produce the luminosity function 
shown in the top panel of Figure [5l Even with the large photometric errors (a ~ 0.07 at 
F814W = 27), the sharp rise of the luminosity function is evident, as is a small contribution 
from brighter stars, which presumably arise from the galaxy's AGB component (more on 
this below). 

To determine the location of the tip of the red giant branch, we applied to our continuous 
luminosity function the edge-detection algorithm 



E{m) = $(m + a{m)) — $(m — a{m)) 



T. 



where a{m) is the photometric error defined via our artificial star experiments (ISakai et al. 



19961 ). This function, which is displayed in the bottom panel of Figure [5l reveals a rather 
wide peak near F814W ~ 27.2. From the figure, the location of the RGB discontinuity is at 
Itrgb = 27.22 ±0.15, where the uncertainty is bas ed on the fu ll-width-half-maximum of the 



distribution. This error is likely conservative fe.g.. ISakai et al.l 19961), though it is similar to 
that expected based on the number of stars available ( Madore fc FreedmanlE995l ). 



To derive the distance to the dwarf, we adopt Mj^trgb = —4.06 ± 0.07 (rand om er- 
ror) as the absolute magnitude of the metal-poor RGB tip (IFerrarese et al.ll2000l ). This 
n umber is very similar to that o f other recent determinations: the empirical calibration 
of ba Costa fc Armandrc^ Jl99oh yields Mf tb/^r = -3.96 ± 0.06 for objects with [Fe/H] = 
—2.3, while a more recent determination by lBellazzini et al.l (120011 ) gives Mj^trgb = — 4.04± 
0.12 for [Fe/H]= —1.7. (The metallicity dependence of Mj^trgb over this range is much 
smaller than these random errors.) Thus we derive (m — M)j = 31.28 ± 0.17 for the 
dwarf. Assuming a foreground re ddening of E(B — V) = 0.025 (ISchlegel et al.l Il998l ) 



(and thus Apsuw ^ Aj = 0.046; ISirianni et al.ll2005l ). we derive a distance modulus of 
(m — M)o = 31.23 ± 0.17, or D = 17. Q ±1.4 Mpc (random error only). This distance is simi- 
lar to the Cepheid, surface brightness fluctuation, and planetary nebula luminosity function 
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distances to the core of Virgo (jJacoby et al.lll990l : iFreedman et al.ll200ll : iTonry et al.ll200ll ). 
Although the three-dimensional structure of the Virgo cluster is complex, the line-of-sight 
depth of the Virgo core (as defined by the early-type galaxies) su rrounding M87 is at least 



2 — 3 Mpc (e.g., iNeilsen fc Tsvetanovl I2OOOI : IJerjen et al.ll2004l ). and the distance to the 



object is consistent with being located near the center of the cluster. 



3.2. Metallicity 



As noted earlier, much of the color spread in the dSph's RGB is due to photometric 
errors. Thus, to measure metallicity, our approach was to derive a mean abundance for the 
galaxy, and then attribute any additional scatter on the RGB to a dispersion in metallicity. 
Since the preceding analysis has shown that the galaxy contains no significant AGB popula- 
tion (also see Section 3.4), we can ignore the alternative that the color spread is due to the 
presence of young or intermediate-age stars. 

To derive the mean metallicity of the dSph galaxy, we compared the ri dge line derive d 
from the CMD in Figure H] with t he scaled-solar abundance isochrones of iGirardil (120051 ) . 
which are the Padova isochrones (iGirardi et al.l I2OOOI . |2002| ) transformed directly to the 
ACS/WFC f il ter sy stem. In addition, we have also used the observed fiducial sequences from 
Brown et al.l (120051 ). which are derived from pho tometry of Mi lky Way clusters covering a 
wide range in metallicity. Figure [6] compares the iGirardil (120051 ) 12.5 Gyr sequence (shifted 
by (m - M)f814W = 31.29 and E(F606W - F814W) = 0.025) to the observed colors of 
the galaxy's RGB stars, both with and without the Am magnitude shifts predicted by 
our artificial star e xperiments. Als o shown for c ompa rison is the M92 fiducial sequence of 
Brown et al.l (120051 ). According to iBrown et al.l (120051 ) . this cluster has an [Fe/H] value of 
—2.14, but is enhanced in a-process elements by [a/Fe] = +0.3, thus implying Z ~ 0.0003. 
Consequently, its position in the color-magnitude diagram is consistent with the Z = 0.0004 
scaled-solar Padova isochrone. 

As Figure E] illustrates, the dSph galaxy is extremely metal-poor: the colors of its RGB 
stars are bluer than the metal-poor M92 fiducial and very near the Z = 0.0001 Padova 
isochrone. For a more quantitative estimate, we can use the mean color of the RGB at 
F814W = 27.6 ±0.1 (or Mf814w = —3.7 at the distance of the dwarf) and interpolate in 
the Padova models; this particular magnitude was chosen to be deep enough to adequately 
sample the RGB (e.g., not at the RGB tip), but not so faint as to be adversely affected by 
large photometric errors and photometric incompleteness. This procedure yields a mean 
metallicity for the system of [M/H] = —2.3 ±0.3, where the quoted error is that of the mean, 
not of the distribution. This error does not include possible systematic errors in the Padova 
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models, but does include a ±0.1 dex uncertainty due to distance. The errors in our color 
calibration (~ 0.03 mag) and the reddening (~ 0.02 mag) contribute an additional ±0.3 dex 
uncertainty to our determination. 

To estimate an upper limit to the metallicity spread in the galaxy, we compared the 
observed width of the RGB to the expected dispersion derived from our artificial star exper- 
iments. To do this, we used the data displayed in the right-hand panel of Figure HJ while 
excluding the one, extremely red (F606W-F814W = 2.06, F814W = 27.5) object from 
consideration. We then formed the sum 

2 _ {Cpbs ~ Cjso) 

i c 

where Cobs is the observed F606W-F814W color of each star with 27.3 < F814W < 28.0, 
Ciso, the predicted color for the star (given its observed F814W magnitude and a system 
metallicity of [M/H] = —2.3), and fXc the expected color dispersion at magnitude F814W, 
derived from our artificial star experiments. The resulting value, x'^/'^ = 1-16 for 95 degrees 
of freedom, indicates that the true red giant branch is slightly wider, but not inconsistent 
with, the distribution expected from a single metallicity population. This same type of 
analysis also demonstrates that any additional (intrinsic) color dispersion on the red giant 
branch must be less than crF606W-F8i4W ~ 0.09 at the 95% confidence level. This limits the 
metallicity spread of our dSph galaxy to < 0.6 dex. 



3.3. AGB Stars and the Age of the Stellar Population 



The presence of stars brighter than the RGB tip have often been used to ascertain 
the existence of intermediate-age (t < 10 Gy r) stars: the more luminous the tip of the 



AGB , the younger the ste llar population (e.g., iMould fc Aaronsonlll982l : lArmandroff et al. 



19931 : ICaldwell et al.lll998l ). Our CMD shows that our dwarf spheroidal galaxy contains very 
few stars brighter than F814W ~ 27.1. However, since the number of stars present in the 
galaxy is sma ll, it is difficult t o make precise measurements of the AGB phase of evolution. 
Moreover, as IWilliams et al.l (l2006al ) have demonstrated, our dwarf galaxy resides within 



a sea of intracluster stars, so any RGB stars foreground to the dwarf can masquerade as 
'brighter' AGB objects. 

Nevertheless, it is possible to place a limit on the presence (or absence) of AGB stars in 
our galaxy. To do this, we used our derived distance to the dSph and the stars' observed colors 
to estimate the bolometric magnitudes of the galaxy's brightest objects. S pecifically, we 
converted the stars' F814W— F606W colors to {V — I) using the color terms of ISirianni et al. 
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( 120051 ). and applied the iDa Costa fc Armandrofj (jl990l ) bolometric correction 

BCi = 0.881 - 0.243(y - /)o (3) 

to generate values of M^o; for all objects fainter than F814W = 26.0. This procedure is 
uncertain at the ~ 0.2 mag level: in addition to ~ 0.05 mag errors associated with the color 
transforms and bolometric correction, our distance to the galaxy is uncertain by ~ 0.15 mag. 
Nevertheless, this calculation is sufficiently accurate for our purpose. 

A histogram of the derived values of Mf,oi for the brightest stars in the right-hand panel 
of Figure H] is shown in Figure [71 Also plotted are the maxim um attainable lumino sities for 
AGB stars produced by 3, 5, 8, and 10 Gyr old populations (iRejkuba et al.ll2006l ). As the 
figure illustrates, our sample of stars contains two objects that may indicate the existence 
of a (small) number of intermediate age objects. Both are sufficiently luminous (~ 0.8 
and ~ 0.4 mag brighter than the RGB tip) to exclude the possibility of image blending. 
However, both are also located on the very outskirts of the the galaxy, where the likelihood 
of contamination is greatest. Since the density of intracluster stars (immediately surrounding 
the dwarf) with 26 < F814W < 27 is 36 arcmin"^, we should expect our sample to contain 
~ 1 contaminating source in this magnitude range. Thus, it is uncertain whether these two 
bright objects are actual members of the galaxy. 

Two of the three remaining stars with —4.2 < Mboi < —4.0 are projected near the 
center of the galax y, and are thus likely to be associated with the dSph. If we apply the 
age-Mhoi relation of lRejkuba et al.l (120061 ) to these objects, then the inferred age of the stellar 
population is more tha n 8-10 Gyr. This is similar to the ages of the old populations seen 
in Local Group dSphs (lMatedll998l ). We cannot exclude the possibility that these objects 
may be blends due to the more crowded environment; if this is indeed the case, this would 
strengthen the conclusion that the galaxy contains no stars with ages less than ~ 10 Gyr. 



3.4. Surface Brightness Profile 



In order to determine the basic structural parameters for our dSph, we performed sur- 
face photometry on smoothed images of our frames. To do this, we began by increasing 
the per-pixel signal-to- noise of the F606W and F814W images by re-binning the data 3x3. 
We then smoothed each frame with a circular Gaussian filter to produce images with a 
~ 0'.'6 point-spread- function, an d fit ellipses to the dSph's c ontours using the el lipse task 
in IRAF/STSDAS (lBuskdll996l . based on the algorithms of IJedrzejewskil (119871 )). To facil- 
itate these fits, the center of the galaxy was kept fixed for each contour, but the position 
angle and ellipticity were allowed to vary to produce the best results. In addition, to aid in 
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the measurement of this low surface brightness object, models were also computed for three 
nearby background galaxies, whose light contaminates the dSph's outer contours. By sub- 
tracting these models from the original image, we were able to improve the fits for the target 
galaxy. Finally, to compute the surface brightness of each fitted contour, a sky background 
was determined using the median pixel value derived within 'empty' regions of our frame. 
The resulting surface brightness profiles (out to a radius of 9'.'3 along the major axis) are 
plotted in Figure [8] as a function of the geometric mean radius (r = ^/ab). 

Figure M confirms that our galaxy is indeed a dwarf spheroidal. The ceritral surface 
brightness of the galaxy is typical of Local Group dwarf spheroidals (IMated Il998l ). with 
fio = 24.54 ± 0.03 in F606W and /xq = 23.70 ± 0.03 in F814W. Mor eover as the dotted line 
in the figure demonstrates, the profile of the galaxy is well fit by a iKingI (119621 ) model. To 
derive this line, we convolved a series of iKing) ( 119621 ) profiles with the smoothed PSF, and 
fit the F814W data using all measurements with errors less than 0.3 mag. This procedure 
yields an excellent fit xV'^ = 7.9/20 = 0.40, with a core radius of = 2'.'6 ± O'.'T and a 
tida l radius of = 10" ± 3" where the uncertainties are derived via Monte Carlo simulations 



see 



Williams et al.ll2006bl ). At the distance of the dwarf, these measurements correspond to 



220 ± 80 pc and rt = 850 ± 250 pc, where the errors include the uncertainties in both 



the fit and the distance. 



We have al so fit our surface brightn ess profiles with the commonly used ISersid (119681 ) 
profile (see also Graham fc Driver 2005), i.e., S = Eoe~^''/^°) '' . Least squares fits to the 
F814W and F606W data yield shape parameters very close to the n ~ 0.5 value expected 
for isothermal distributions: for the F814W image. Eg = 23.70 ± 0.03 mag arcsec"^, tq = 
3'.'16±0'.'16 and n = 0.61 ±0.12, while the F606W profile has Sq = 24.52±0.03 mag arcsec^^ 
ro = 3'.'34 ± 0'.'06 and n = 0.59 ± 0.12. Both fits (also plotted in Figure [HD are only slightly 
worse than those for the King model [yd 1'^ = 0-59 and 0.63 for the F814W and F606W fits, 
respectively), but are still a very good match to the data. 

We can also use the data of Figure [8] to compute the integrated magnitude of our galaxy. 
If we sum the flux contained within each elliptical isophote, then the total magnitude of the 
dwarf spheroidal is F606Wtot = 20.56 ± 0.05 and F814Wtot = 19.86 ± 0.05. These values 
are consistent with the F606W = 20.46 and F814W = 19.74 magnitudes obtained from 
simple photometry using an 11" circular aperture centered on the galaxy. T o conv ert these 
magnitudes to V and /, we used the color transformations of ISirianni et al.l (120051 ) . and to 
obtain absolute magnitudes, we applied the distance and reddening values given in Section 
3.1; no additional corrections to our integrated magnitudes were required as our profiles 
extend to very faint surface brightnesses. The data imply that the dSph galaxy has a total 
absolute magnitude of My = —10.6 ± 0.2 (and Mj = —11.4 ± 0.2), where the error again 
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includes both the uncertainties of the photometry and the distance. 



3.5. Comparison with Local Group dSph Galaxies 



Detailed observations in the Local Group demonstrat e that dSph galaxies possess a 
varie ty of star-formation histories and chemical abundances ( iGallagher fc Wyselll994l : [Mateo 
19981 ). This has been attributed to the effect of environment on their formation and evolution. 
If so, then the properties of such galaxies in dense clusters might be distinctly different from 
those of dSphs seen locally. Unfortunately, althou gh a few dwarf galax i es have been studied 



in other nearby group s, such as the M81 system (jCaldwell et al.lll998l : iKarachenstsev et al. 



20051 : iDa Costal l2005l ). very little is known about dSphs in rich clusters. 



Figure [9] compares the observed central surface brightness, core radius, and metallicity 
of our Virgo dSph with the corresponding properties of dSph galaxies in the Local and 
M81 Groups. The most striking aspect of the figure is how normal our dSph appears. The 
galaxy's structural parameters are similar to those of the Leo II and Ursa Minor dwarfs, i.e., 
ty pical of objects present in the l ocal n eighb orhood. This res ult is consistent with the studies 
of iDurrelll (119971 ) , iHilker et al.l (120031 ) , and ICaldwelll (120061 ) , who also found no significant 
structural differences between the dwar fs of the Local G roup and those of Virgo and Fornax. 
For comparison, the Virgo results from ICaldwelll (120061 ) are also plotted in Figure O 



We reach a similar conclusion from the metallicity measurement of our dSph galaxy. 
Dwarf elliptical and spher oidal galaxies in the Local Group obey a well-established luminosity 



metalhcity relation (e.g., ISmithlll985l : ICaldwell et al.lll992l . Il998l : ICote et al.ll2000l : ICaldwell 



20061 ). where the more luminous galaxies are (on average) more metal-rich. This result, 
which is usually at tributed to increased mass loss in systems with small potential wells (e.g., 
Dekel fc Silklll986l ). predicts the metallicity of our dSph galaxy fairly accurately. However, 
galaxies in clusters are expected to lose much of their mass to int racluster space via tidal 



encounters with o ther galaxies and with the cluster as a whole (e.g.. lMoore et al.lll996l . Il998 
Bekki et al.ll200ll ). If our dSph had been exposed to these forces for any length of time, 
we would expect its metallicity to be more in line with that of a higher-luminosity object. 
This does not appear the case: if anything, our galaxy's metallicity lies below the mean 
luminosity-metallicity relation. This suggests that most of the galaxy's original stars are 
still bound to the system, or perhaps the luminosity-metallicity relation differs in the cluster 
environment. 



The structure of the galaxy appears to support this conclusi on. T i dal fe atures are of- 
ten detected in dSph galaxies via photometric deviations from a iKingj (119621 ) profile (e.g.. 
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Irwin fc Hatzidimitrioij Il995l : iPalma et al.l l2003l : iMajewski et al.l |2005| ) . As Figure [8] demon- 
strates, our Virgo dwarf is well-fit by a King model all the way down to fiy ~ 27.5 
[fij ~ 26.5), where the galaxy's surface brightness merges into that of the intracluster 
background (jWilliams et al.l l2006al ). Moreover, a plot of {x,y) positions using only those 
stars with colors and magnitudes similar to those of the galaxy's RGB objects (i.e., 0.4 < 
F606W - F814W < 1.3 and 26.8 < F814W < 28.0), reveals no obvious evidence for tidal 
distortion (Figure [T0|) . We note that our observations are not deep enough to see any low- 
level tidal disruption of this dSph galaxy: the large streams see n in the Coma and Centaurus 
Clusters are only vis i ble at surface brightnesses bel ow nv ~ 26 (ITrentham fc Mobasherlll998 



Gregg fc WestI Il998l : ICalcaneo-Roldan et al.l l2000l ) , and the tails associated with the Local 



Grou p galaxies Carina and Ursa Minor are at /ly > 29 (jPalma et al.ll2003l : iMajewski et al. 



20051 ). Thus we cannot unambiguously state whether the dSph galaxy has undergone any 
tidal stripping. Nevertheless, the metallicity arguments above suggest that the observed 
dSph is not the low-luminosity remains of a much more massive system (e.g., as expected in 
galaxy harassment scenarios). 

Tidal threshing is not the only process that could affect dwarf galaxies in the cluster en- 
vironment. In th e Local Group, rarn pressure stripping rnay be a prime factor in the evolution 



of dS ph galaxies (iLin &: Faberlll983l : Ivan den Berghlll994j : Grebel et al.ll2003l : iMarcolini et al. 
2003h . Indeed, the N-body simulations of iMaver et"airj2006h suggest that both tides and 



ram pressure are needed to explain the properties of Local Group objects. If true, then the 
effects of intracluster gas should be magnified in the denser regions of the Virgo Cluster core. 
If low- luminosity dwarfs can barely hold onto their gas (and enrich themselves) in the Local 
Group, then Virgo dSphs should easily be stripped by the Virgo intracluster medium. The 
result would be a dSph galaxy much like those observed in the Local Group (as suggested by 
the structural similarities), but with less chemical enrichment, and a very narrow red giant 
branch. 

Our dSph does show some evidence for enhanced pressure stripping: as Figure [9] shows, 
the galaxy's mean metallicity is slightly less than that predicted from dSph observations in 
the Local and M81 Groups. However, our RGB photometry is not precise enough to test 
for self-enrich ment. Local Group dwarfs typically have metallicity dispersions of > 0.3 dex 
( lMatedll998l ). Our data, though consistent with a zero abundance spread, is best fit with 
an RGB color dispersion of ctf606W-F814W = 0.05 (~ 0.4 dex) and still admits a metallicity 
spread as great as crF606W-F8i4W = 0.09 (~ 0.6 dex). Thus, while the data are suggestive, we 
cannot say for certain that the Virgo environment has had any additional effect on the dSph's 
chemical properties. Deeper photometric studies of other Virgo dSph galaxies would be very 
useful in ascertaining any differences between the chemical evolution of low-luminosity dwarf 
galaxies in clusters and those in the group environment. 
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Discussion 



From the analyses in the previous section, it is clear that the dSph galaxy in our field 
lies in the Virgo cluster and appears relatively undisturbed by its surroundings. Its stellar 
population is very metal-poor ([M/H] ~ —2.3) and old, with no indication of stars younger 
than ~ 10 Gyr. In this sense it is similar to the purely old Local Group dSph g alaxies of 
Ursa M inor and Draco. It is also similar to the dwarf galaxies recently observed by ICaldwell 
( I2OO6I ) in the Virgo Cluster core - their central surface brightnesses and metallicities are also 
included in Figure [H They, too, obey the luminosity-metallicity relation for small galaxies 
and appear undisturbed by the cluster environment. 

The appearance of such 'pristine' old, low-luminosity dwarf galaxies within the Virgo 
Cluster begs a question - how do these objects last so long in the cluster? S mall galaxies 



should be easy prey for the tidal forces associated with the cluster environment (IMoore et al. 



19981 ). and, indeed, the velocity distribution of Virgo dE galaxies supports the idea that at 
least some fraction of the cl uster's dwarfs are the stripped remains of larger, late-type systems 
(e.g., IConselice et"aDl200lh. Furthermo re, it is plausible that many of Virgo's metal-poor 
intracluster stars (jWilliams et al.ll2006al ) originated in such galaxies. Nevertheless, it seems 
clear that some dwarfs remain intact; perhaps these are dwarfs that have recently fallen into 
the cluster, or have rather high M/L ratios. 

A weak constr aint on the history of our dSph galaxy can be obtained from its observed 
tidal radius. From iKingj (119621 ). the limiting radius of the galaxy should be 



Tt = Rp 



Mc(3 + e) 



1/3 



(4) 



where Mg is the mass of the galaxy, Mc, the mass of the cluster, Rp, the galaxy's pericluster 
radius, and e, the eccentricity of the orbit . The mass o f the galaxy can be obtained by analogy 
with Local Group objects: according to (lMatedll998l ). dwarf galaxies with My ~ —10.6 have 
mass-to-light ratios of M/Ly ~ 10, implying Mg ~ 1.5 x lO'^Mp,. The V ir go clu ster mass is 
similarly obtainable from the x-ray observations of iNulsen fc Bohringerl (119951 ): their data 



give a core radius of 56 kpc, a mass per unit length of 12 x lO^^M© kpc~^, and a total mass 
within 180 kpc (the projected cluster distance of the galaxy) of 2 x lO^^M©. A value of 
e ~ 0.5 then excludes any orbit which extends into the inner ~ 80 kpc of the cluster. 

Of course, there are any number of ways to satisfy this condition. For instance, the 
dSph may just now be falling into Virgo for the first time. Cepheids, the planetary nebula 
luminosity function, the surface brightness fiuctuation method, and the tip of the red giant 



branch, all place the core of Virgo at a distance of betweeri 14.5 and 17 Mpc (IFreedman et al. 



2001; 


Jacobv et al. 


1990; 


Tonrv et al. 


2001; 


Jensen et al. 


2003 



2004 



Caldwell 
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20061 ). Our dSph distance estimate of 17.6 ± 1.4 Mpc is consistent with this range of values, 
but it also admits the possibility that the galaxy is slightly behind the main body of the 
cluster. If so, then the dSph's location may associate it with the infalling M86 system, 



which is ~ 3 


Mpc 


Jerjen et al. 


2004) 



smaller cluster /group, and the shredding of the galaxy is only now about to begin. 

If the dwarf did originate and live much of its life in the cluster environment, then it 
would have to have a high mass to withstand the tidal forces of Virgo's densest regions. 
This would require a mass-to-light ratio that is much larger than the 'normal' value of 
M/Ly ~ 10 observed in the Local Group. CDM models do suggest that dwarf galaxies can 
form out of dark matter ha los much more ma s sive than dynamical measurements sugges t 
(i.e., M ~ fe.g.. Istoehr et allbooi Iflavashi et aPboOsI : Ecravtsov et al.ll2004[ ). 



and such objects would be able to resist tidal disruption. However, galaxies with large dark 
matter halos would no t easily be stripped of their ISM by ram pressure or galactic winds. 
Mori fc BurkertI (120001 ) point out that while low mass (M ~ IO^-'^Mq) objects will rapidly 
lose their gas to the intracluster medium, more massive systems are likely to hold on to 
at least some of their interstellar medium. As a result, their stellar populations should be 
younger and more metal-rich than their lower-mass counterparts. While it is beyond the 
scope of this paper to investigate the precise masses and conditions under which early gas 
removal is expected, the old, metal-poor population of our dwarf spheroidal argues against 
the presence of a very massive dark halo. The existence of such a massive halo could, 
however, be tested via high-resolution spectroscopic observations with the next generation 
of large telescopes. 

Unfortunately, our lack of knowledge about the dwarf's current location and orbital 
characteristics precludes our placing a strong constraint on its origin. While we have no 
data on the presence (or absence) of H I gas in the galaxy, that gas has certainly not created 
many new stars in the past ~ 8 Gyr. Thus, if gas were detected, it would not affect our 
basic conclusions. 



5. Summary 

We report the discovery of a dwarf spheroidal galaxy (with My = —10.6 ± 0.2) on deep 
ACS F606W and F814W images of an intracluster field in the Virgo Cluster. The distance to 
the galaxy is 17.6 ± 1.4 Mpc based on the location of its TRGB; this value is consistent with 
a location close to the core of the Virgo cluster, although we cannot rule out the possibility 
that this object is part of the M86 subcluster falling into Virgo from behind. The galaxy is 
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clearly resolved into stars, and our observations extend more than a magnitude down its red 
giant branch. 

We find that the galaxy is composed entirely of old (> 8 — 10 Gyr), and is very metal- 
poor, with [M/H]= —2.3 ±0.3. This metallicity is consistent with that expected for a galaxy 
of its luminosity; if anything, the galaxy lies slightly below the mean luminosity-metallicity 
relation. Thus, the object is not the remains of a larger galaxy that has been tidally stripped 
or harassed. Moreover, the dSph's structural properties are similar to those derived for dwarf 
galaxies in the Local and M81 groups, suggesting that many of the same physical processes 
that govern the formation and evolution of local dSphs also apply to this Virgo Cluster 
object. Based on this similarity, we suggest that our dSph galaxy is likely on its initial infall 
into the center of the cluster. 
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Fig. 1. — Color image of a 28" x 24" section of our image, centered on the dSph galaxy. In 
the image, blue represents 2x F606W— F814W, green represents F606W, and red represents 
F814W. The ellipse denotes the boundary used to defined a subsample of stars which 
minimizes contamination (see text). North is to the top, and East is to the left. 
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Fig. 2. — Fraction / of recovered artificial stars for the F814W (filled circles) and the 
F606W (open circles) images. The F606W data is based on added stars with input col- 
ors F606W-F814W=2.0, and that of the F814W data is from stars with input colors 
F606W-F814W=0.0. 
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Fig. 3. — Output color-magnitude diagram for our artificial stars. The dashed line denotes 
the input colors for the stars (F606W— F814W) = 1.00. Note the slight offset in the colors at 
faint magnitudes. This is due to the different limiting magnitudes of the two filters. 
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Fig. 4. — (a) The color-magnitude diagram (in the VEGAMAG system) for the 611 stellar 
objects located in a 66" x 48" region centered on the dSph galaxy, (b) The 'dwarf-only' 
CMD, formed from a subset of 181 stars located within the inner elliptical region shown in 
Figure [H The dotted lines denote the 50% completeness levels, while the error bars represent 
the typical photometric uncertainties. Note the discontinuity at F814W ~ 27.1; this is the 
tip of the red giant branch. 
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Fig. 5. — The top panel shows a continuous luminosity function formed by co-adding the 
Gaussian representation of stars in our minimally contaminated subsample (see text). Ob- 
jects with F606W— F814W > 1.3 have been exchided from the analysis. The bottom panel 
is the result of applying the Sobel edge-detection algorithm. This panel suggests that the 
tip of the red giant branch is at FSUWtrgb = 27.22 ± 0.15. 
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Fig. 6. — A comparison of observed red giant branch of our dSph galaxy with the 12.5 Gyr 
models of lCirardil J2005h . The solids are (from left to right) Z = 0.0001, 004, and 0.001. 
The dashed hnes show the M92 fiducial (Z ~ 0.0003) from ISrown et all J2005h . All the 
curves have been shifted to a distance modulus of (m — M)fsi4W = 31.29 and reddened by 
£'(F606W — F814W) = 0.025. The points represent the observed mean red giant branch, 
while the error bars show the uncertainties in the mean. The left panel displays the raw 
measurements; the right panel shifts the observed colors by an amount predicted by our 
artificial star experiments. 
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Fig. 7. — A histogram of bolometric magnitudes for the brightest stars in our minimally 
contaminated dSph subsample; no corrections for photometric incompleteness have been 
applied. The arrows denote the maximum attainable brightnesses of AGB stars forme d from 
populations with ages of 3, 5, 8, and 10 Gyr (from Figure 19 of iRejkuba et al.ll2006l ). The 
two brightest stars are likely interlopers. The data suggest that the dSph population is at 
least 8 Gyr old. 
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Fig. 8. — F814W and F606W surface brightness profiles derived by fitting elliptical contours 
to the smoothed images of the dSph galaxy. T he data are plotted as a function of the 
geometric mean radius r = ^/ab. The best-fitting Kinej ( 
and Tt = 10" ± 3" is shown as a solid line. The iKingl (j 
X^/z^ = 0.40. The dashed lines denote the best fits to a ISersid (119681 ) profile; see text for 
details. 
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Fig. 9. — A comparison of the structural and chemical properties of the Virgo dSph (filled cir- 
cle) with measurements for the dwarf spheroidals in the Local and M81 Groups (open circles). 
The top panel shows central surface brightness, the middle panel core radius, and the bot- 
tom panel metallicity, with the error bars representing the metallicity dispersion. The filled 
triang les denote measurements of faint dE/dSph galaxies in t he Virgo Cluster fromlCaldwell 
(120061). The Local G roup data come from the compilati ons of llrwin Sz Hatzidimitrioul (Il995l ). 



Grebel et al.l f 20031 ) and McConnachie fc IrwinI J2006h. with additional information for in- 



dividual dwarfs from Saviane et al. 



1996h.|Palma et al.l J2003h . and iHarbeck et all J2005h . 



The M81 dwarf data come from lCaldwell et al.l (119981 ). Note that the properties of the Virgo 
dwarf are completely consistent with those of dwarfs in these low-density systems. 
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Fig. 10. — Spatial distribution (in pixels) for those stellar objects near our dwarf galaxy with 
the colors and magnitudes of metal-poor red giant stars (i.e., 0.4 < F606W — F814W < 1.3 
and 26.8 < F814W < 28.0). Each pixels represents 0'/03, or ~ 2.2 pc at the distance of 
Virgo. There is no obvious evidence for the tidal shredding of the dwarf galaxy. 



